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Abstract. Intercalation of polyatomic molecules into a superconductor can drastically affect the properties of the com- 
pound. A mechanism responsible for a large increase in T c for such systems is proposed. It explains the recent remark- 
able observation of high T c superconductivity in the hole-doped C60/CHX3 (X=Cl,Br) compounds and the large shift 
in their T c upon CI— >Br substitution. The increase in T c is due to contribution to the pairing arising from the interaction 
of electrons with the vibrational manifold of the molecule. The proposed mechanism opens up the possibility to observe 
a site-selective isotope effect. We also suggest that intercalating CHI3 would further increase the critical temperature to 
T C ~140K. 

PACS. 74.70.Wz Superconductivity; Fullerenes and related materials - 74.72.-h High-Tc compounds 



This paper is concerned with the impact of intercalation by 
polyatomic molecules on the properties of superconductors. 
The study has been motivated by the recent remarkable ob- 
servation of high T c caused by intercalation of CHBr3 (bro- 
moform) and CHCI3 (chloroform) molecules into hole-doped 
fullerides [fij (see also description of [Jl|] in Ref. [^]). The dis- 
covery of high temperature superconductivity [|lj| in these sys- 
tems (T c = 1 17K for C 60 /CHBr 3 ; T c = 80K for C 60 /CHC1 3 ) has 
attracted a lot of attention and raises the fundamental question 
about the nature of this phenomenon. 

The usual superconducting fullerides (see ,e.g., the reviews [Q, 
Q]) are chemically electron-doped compounds (e.g., Rb3C6o, 
KsCgo)- It is believed that their superconductivity is due to the 
coupling of electrons to the intramolecular vibrational modes. 
The vibrational spectrum spreads over a broad region (Q.l ~ 
270cm~ 1 — * Q.h — 1500cm -1 ). Strictly speaking, all modes 
contribute to the pairing but the question of which region of 
the vibrational spectrum plays a major role is still controversial 
(see, e.g. @]). 

A new exciting development was described in [g]; with the use 
of gate-induced doping, a hole-doped fulleride compound was 
created without chemical dopant. As a result a major increase 
in T c to 52K(!) was observed. This is caused by the different 
properties of the valence band relative to the conduction band 
such as the density of states. Recently, this technique was com- 
plemented by intercalating the molecules CHBr3 and CHCI3 
into the crystal [[j]]. These molecules are placed in regions near 
the C60 units. This has resulted in a dramatically higher T c 
again: T c ~ 117K(!). It was indicated in [|J that the change in 
lattice spacing may be responsible for the increase in T c . This is 
an interesting scenario which needs to be studied in detail. Here 
we propose another additional mechanism. It is of interest in its 
own right; furthermore, it can provide for a quantitative expla- 
nation of the drastic increase in T c observed in [fit]. 



It is very essential to note that the experimentalists [Uj] observed 
not only a remarkable high T c for the C60/CHX3 compound 
(X= Br,Cl), but also a large difference in the value of T c for 
X= Br vs. X= CI. Whereas T c = T^ r ~ 117K was observed 
for X=Br, the value of T c for X= CI turned out to be much 
smaller (~ 80K), although still very high. The structures of the 
bromoform and chloroform molecules are similar and therefore 
a microscopic explanation of the strong change in the critical 
temperature is also of definite interest. 

The picture we propose is that the intercalated molecules 
themselves actively, and importantly, participate in the pairing. 
The internal vibrational modes of the molecules provide for an 
additional attraction which leads to higher T c , Such a mech- 
anism was considered by one of the authors in [Q], Therefore, 
our major focus is on the vibrational spectra of these molecules. 

Let us start with the parent hole-doped fulleride (T c = 52K). 
The order parameter for this superconductor is described by the 
equations: 



A(co„)Z(co„) 



Jir£[XiZ)((D B -«V; 61 )-/;*] 



A(co„/ 



Z(co„) = l+A,i7tr^D(co„-co„/;ni) 



C0„/ 



(1) 



(2) 



Here X\ is the electron-phonon coupling constant, Cli is the 
characteric (average) intramolecular vibration frequency 

(fil =< i2 2 see, e.g., Ref. [|), D = £l\ [Cl] + (co„ - ov) 2 ] 

is the phonon Green's function, co„ = (2« + 1 )kT, T = T c , Z is 
the renormalization function, and // is the Coulomb pseudopo- 
tential. We are employing the temperature Green's functions 
formalism (see, e.g 
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Eqs. ([j]), (g) correlate the value of T c with the values of the 
parameters X\, fi*. It is natural to assume that hole-doped 
fullerides differ from electron-doped compounds obtained 
by chemical doping in the magnitude of the coupling constant 
X\, whereas the values of the characteristic phonon frequency 
Cli and /j* should remain similar. Following Ref. [0], let us take 
as an estimate Cl\ ~ 0.1 eV and /j* = 0.15. Then, the value of T c 
= 52K observed in [||] implies the value X\ w 0.75. We would 
like to emphasize that if we were to adopt a different value 
of Cli (e.g., £2i = 80meV, or even 50meV), this would affect 
some intermediate parameters (see below), but the final result, 
namely the value of T c for the C6o/CHBr3 compound, would 
not be affected in any noticeable way. 

In the intercalated Cgo compound, polyatomic molecules 
CHX3 (X=Cl,Br) are added at interstitial sites, (see, e.g., Fig. 1 
in Ref. [Q]). As indicated above, their presence is found to result 
in two phenomena: 1. A remarkable high T c (T c ~ 117K) for 
X=Br; 2. A very high, but significantly lower T c (~ 80K) for 
intercalation by chloroform (X=C1). 

The molecules have a relatively complex vibrational spec- 
tra (see, e.g., Ref. [j^]), but it can be modeled as an additional 
phonon mode. Then, instead of equations (fil ^ we have: 



A(co„)Z(co„) = 



1=1,2 



A(co„ 



CO,,/ 



CO,/ 



z(co„) = i+7tr V itYDinn-n^Qi)^ 



(3) 



(4) 



The values of X\ and £2i correspond to the parent, non-interca- 
lated state [see Eq. (0)]. We assume that the value of /j* is also 
essentially unchanged, although the presence of molecules and 
their electrons might lead to an effective decrease in /u* (see 
below). 

The presence of the second term (~ X2) m the sum on the 
right-hand side of Eq. (Q) leads to an increase in T c relative to 
that of the parent compound (T° ~ 52K). The magnitude of this 
increase depends on the values of £^2 and Xi- 

As is known, the value of the coupling constant X is de- 
scribed by the McMillan relation [Pj: 



X 



v<I z > 
M < Q? > ' 



(5) 



where v is the density of states, / is the matrix element, and M 
and £2 describe the vibrational spectrum. We will show that the 
differences in the vibrational spectra of the CHX3 molecules 
can account quantitatively for the observed large difference in 
T c upon chemical substitution (Cl^Br). Let us now turn to 
quantitative analysis. 

The vibrational spectrum of the CHCI3 molecule (see, e.g. 
||8[) contains five major modes which can contribute to the pair- 
ing. These modes are Q.^ cl ~ 262cm -1 , D.^ cl ~ 366cm -1 , 
q(3)cz ^ 668cm- 1 , £2 (4 )° ~ 761cm- 1 , £2< 5 ) c/ ~ 1216cm" 1 
(the notation £2 C/ means that we consider the chloroform mole- 
cule; the modes Q.W are degenerate). All these 
modesjire Raman active and can couple to the conduction elec- 
trons 



Table 1. Calculation of Xf, Xf and T c for C 60 /CHBr 3 . X\ : coupling 
to intramolecular vibrations of C60, X^: additional coupling to the vi- 
brational modes of CHCI3 or CHBr 3 . The coupling constant X2 for 
CHBr 3 is fixed by the other set of parameters (see text). 





hole-doped C50 


C 60 /CHC1 3 


C 60 /CHBr 3 


Xi 


0.75 


0.75 


0.75 


h 





0.2 


0.55 


T c [K] 


52" 


80 


110 



Ref. 5h 



Consider now the CHBr3 molecule. It has the same configu- 
ration and a similar set of modes, but all of them are noticeably 
softened. Indeed, according to Ref. [§], Q.W Br ~ 154cm -1 , 
a (2)Br _ 222cm -1 , Q.^ Br ~ 539cm -1 , Q.W Br ~ 656cm -1 , 
Q.( 5 > Br ~ 1142cm -1 . As a result, it can be expected that the 
well-known "softening" mechanism (cf., e.g. [|l0||) manifests 
itself. That is, T c is expected to increase as a result of the sub- 
stitution. 

The effective coupling constants [see Eq. (5)] also contain 
the mass (M) factor. It is crucial that all these modes in CHX3 
molecule correspond, mainly, to the motion of the light car- 
bon ion relative to the heavy X3 (X= CI, Br) unit. As a result, 
M w Mq and the difference in coupling strength between the 
two compounds is mainly determined by the aforementioned 
frequency shift (Xi ~ n ; -2 ). Note that the hydrogen stretching 
mode is not essential since it lies much higher in energy and, if 
desired, can be singled out by isotope substitution. 

Let us consider the C60/CHCI3 system first. Usingthe above 
mentioned vibrational modes taken from literature Oq], we de- 
termine the average frequency to be £^2 = Qj 1 = 85meV (care 
should be taken of the degeneracy of the modes). Performing 
a numerical calculation with Eqs. (§),(§) and the derived fre- 
quency £22, we obtain that the observed value T c = 80K [[[]] im- 
plies the value X^ 1 ~ 0.2. Such a value is small and perfectly re- 
alistic, confirming that the observed rise upon intercalation can 
be assigned to the action of the added intramolecular modes. 

Let us now substitute Cl^Br and compare C60/CHCI3 and 
C6o/CHBr3. This comparison is an important step in our anal- 
ysis. As was stressed above, the total coupling constant X Br 
changes relative to X% 1 primarily as a result of the vibrational 
frequency shifts. Using the literature values for the vibrational 
frequencies Q.f r and flf 1 (see above), we calculate that X Br ~ 
2.65A-2'. Therefore, the "softening", indeed, leads to a notice- 
able increase in the coupling constant. 

0.55. Using 



Since X^ 1 



■ 0.2 (see above), we obtain A-f' 



this value for X Br and the frequency £2f r — 70me V infered from 
the above vibrational spectra, we solve Eqs. (Q) and (Q) and 
obtain T c ~ IIQK for the C6o/CHBr3 compound. It is impor- 
tant to emphasize once more that during this last step the value 
of T Br has been obtained without any additional adjustable pa- 
rameter. Table [j] summarizes our numerical results. The derived 
value T c ~ 11 OK is close to the experimentally observed one 
T c ~ 117K [JTJ . Since our calculation made use of some ap- 
proximations, the agreement can be considered as rather good. 

Refs. [|l[| and [|TJ also describe the study of electron-doped 
systems. We first note that doped electrons are filled into the 
conduction band made of orbitals with T\ u symmetry (at the 
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r point), whereas doped holes fill in the valence band which 
has H\ u symmetry (see, e.g., Ref. [Q.) A group theoretical 
analysis shows that in the former case, electrons couple to in- 
tramolecular vibrations of the C60 cluster belonging to the five- 
fold degenerate h g irreducible representation. In the hole-doped 
case, on the other hand, the higher symmetry of the valence 
band allows for coupling of holes to additional vibrational modes. 
Both h g and the four-fold g g vibrations are Jahn-Teller active in 
hole-doped systems. One can therefore expect that the overall 
coupling is larger for holes than for electrons, thereby enhanc- 
ing T c . 

The procedure described above for hole-doped C60/CHX3 
can be applied to the electron-doped case and reproduce the 
results obtained in Ref. [[j]] as well. Indeed, from Eqs. (|l]j^) 
we conclude that the observed value T c = 1 IK Jl~T] ] is reached 
for A-i w 0.5. Since the frequencies of the h g and g g modes 
are in the same range Jl3[], we take the same value of Cl\ for 
electron- as for hole-doped systems. From Eqs. we further 
obtain the experimentally observed T c = 18K for C60/CHCI3 
with an additional coupling X^ 1 ~ 0.07. The scaling relation 
Xf r ~ 2.65A-2' then leads to A.f r ~ 0.2. As a result we obtain 
T c w 30K for electron-doped C6o/CHBr3 . This is in good agree- 
ment with the experimental value T c ~ 26K ^ . 

We have focused on the vibrational spectra of the inter- 
calated molecules. Note that an additional contribution to the 
pairing could also come from the Coulomb interaction and the 
corresponding polarization of the electronic systems of the mo- 
lecules, virtual transitions between electronic levels (see Ref. 
1 15]). This additional attraction would correspond to the elec- 
tronic energy scale and can be described as an effective de- 
crease of /u* (cf. Ref. ||l(|). Note that the polarizability a Br = 
(XcHBr 3 is much larger than a c/ Jl7|]. This question will be con- 
sidered elsewhere. Another factor indicated in [[!]] is the in- 
crease in "ionic" radius upon intercalation and Cl^Br substitu- 
tion. However, recent experiments indicate that lattice stretch- 
ing resulting from the intercalation cannot account for the in- 
crease in T c [[Til], It thus remains an open question whether this 
effect by itself is capable of producing such a giant impact on 
T c . 

We propose also the following experiment directly related 
to our approach. Namely, one can observe site selective isotope 
effect by substituting 12 C^ 13 C on the intercalated molecules. 
Such substitution will affect T c because the vibrational mani- 
fold of these molecules contribute to the pairing. From Eqs. (|3jQ) 
and the experimentally observed shift in D.^ cl j = 1 , . . . , 5 for 
CHCI3 we obtain that AT C = T c ( l2 C) - T C ( 13 C) ~ -0.2K. 
Though small, such a shift can be measured experimentally 
(cf. [|o[). The corresponding isotope coefficient 
a£' = -(M C /AM C )(AT C /T C ) is of order 0.03. If, in addition, 
one also substitutes 35 C1^ 37 C1 a shift in T c of about ~ — 0.35K 
can be observed. 

A more precise numerical treatment would require proper 
use of tunneling spectroscopy, similar to that described in Ref. 
[14], in order to obtain the accurate form of <x 2 (Q.)F(Q.) 0. Our 
aim was to identify the mechanism responsible for the strong 



1 Note that the analysis of a 2 F has to be made with considerable 
care, as the vibrational manifold of the intercalated molecules partially 
overlaps with that of the Ceo intramolecular modes. 



increase in T c for the intercalated compounds. We argued that 
this mechanism involves the vibrational manifolds of the inter- 
calated molecules. We demonstrated that this mechanism can 
naturally lead to high values of T c , as observed in Ref. [[!]]. 

Finally, we would like to propose to carry out the same 
experiment as in Ref. [[j]] with the compound C60/CHI3. The 
corresponding modes are Q^ 1 ~ 1 10cm _1 , Q^ 1 ~ 154cm- 1 , 
q(3)/ ~ 425cm- 1 , Q( 4 ) ; ~ 578cm- 1 , Q.^' ~ 1068cm- 1 [0. 
Our analysis shows that such intercalation leads to T c ~ 140K(!) 
Note also that the site selective isotope effect 12 C<-> 13 C for 
CHI3 molecules has the additional advantage that iodine has 
no natural isotopes. 

The mechanism of pairing proposed here and connected 
with the interaction of electrons with vibrational manifolds of 
intercalated molecules can provide for a drastic increase in T c . 
The search for new systems which would similarly permit the 
intercalation by appropriate polyatomic molecules is an excit- 
ing and promising direction in the quest of new high tempera- 
ture superconductors. 

We thank B. Batlogg for sending the manuscript (Ref. [jl|]) prior to 
publication. 



Note added in proof: 

Recent tunneling measurements [^] have shown the appear- 
ance of additional peaks upon intercalation. This additional 
structure corresponds to the vibrational manifold of the molecules, 
and Cl^Br substitution leads to the "softening". The data of 
Ref. p2y support our theory. 
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